Adiabatic cooling of a tunable Bose-Fermi mixture in an optical lattice 



O. S0e S0rensen and N. Nygaard 
Lundbeck Foundation Theoretical Center for Quantum System Research, 
Department of Physics and Astronomy, University of Aarhus, DK-8000 Arhus C, Denmark 



o 
o 

(N 



o 



I 

O 

o 



> 

in 

rn 

O 
On 
O 



X 



p. B. Blakie 

Jack Dodd Centre for Photonics and Ultra-Cold Atoms, 
Department of Physics, University of Otago, New Zealand 

We consider an atomic Fermi gas confined in a uniform optical lattice potential, where the atoms 
can pair into molecules via a magnetic field controlled narrow Feshbach resonance. Thus by adjusting 
the magnetic field the portion of fermionic and bosonic particles in the system can be continuously 
varied. We analyze the statistical mechanics of this system and consider the interplay of the lattice 
physics with the atom-molecule conversion. We study the entropic behavior of the system and char- 
acterize the temperature changes that occur during adiabatic ramps across the Feshbach resonance. 
We show that an appropriate choice of filling fraction can be used to reduce the system temperature 
during such ramps. 



I. INTRODUCTION 

One of the most fascinating consequences of many- 
body quantum theory is the striking difference in be- 
havior between systems made from identical bosons or 
fermions at low temperatures. The use of Feshbach reso- 
nances in atomic systems is now a routine technique for 
pairing atoms into diatomic molecules [J, 0, 0, 3 S, 0, 
0,01, and leads to the intriguing scenario: the controlled 
conversion of an ultra-cold degenerate Fermi gas into a 
system of bosonic dimers. Hence, in experiments it is 
possible to control the statistics of the system by tuning 
across the Feshbach resonance. In this paper we consider 
a particular case of such a system confined in an optical 
lattice, schematically shown in Figure [H Atom- molecule 
conversion in an optical lattice has been demonstrated 
in several recent experiments [l^, [ill, [l3, [l3, [l3| • We 
develop formalism to quantify the effect that ramping 
across the Feshbach resonance has on the temperature of 
the system in the adiabatic limit by developing expres- 
sions for the entropy in the nearly degenerate regime. 

The essence of a Feshbach resonance is that free atom 
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Figure 1: Schematic diagram of the system under consider- 
ation: Fermionic atoms in two different spin states occupy 
states in an optical lattice and are coupled by a Feshbach 
resonance into bosonic dimers. Because the bosonic dimers 
have a larger polarizability they experience a deeper lattice 
potential with a different spectrum to that for the atoms. 



pairs are coupled to a discrete dimer state in a closed 
scattering channel. The energy of this closed channel 
bound state is detuned from the atomic threshold by an 
amount i?rcs, which is tunable with an applied magnetic 
field, since the unbound atoms and the closed channel 
dimer have different magnetic moments (e.g. see [lB|)- If 
ii^ics < the closed channel bound state corresponds to 
a bound molecular state of the coupled system. When 
Eres > the discrete state is embedded in the atom- 
ic continuum and has a finite lifetime. Consequently, it 
manifests itself not as a true bound state, but as a scat- 
tering resonance. However, if the width of the Feshbach 
resonance is sufficiently narrow, the lifetime of the closed 
channel dimer state is long enough that we may consider 
it a quasi-bound state of the system. Hence, character- 
izing both the bound and the quasi-bound state of an 
atom pair as molecules, the number of unbound atoms 
and molecules in the gas are well-defined at any given 
instant. 

The thermodynamics of such a Feshbach-resonant Fer- 
mi gas have been studied in free space [lB| as well as in 
a harmonic confinement (itI . [l^ . and in previous work 
we have considered the case of a lattice potential and 
analyzed the behaviour of the chemical potential and 
the atomic and molecular populations in the degener- 
ate regime (/cbT < O.lEp). There has also been a great 
deal of work on systems with broad Feshbach resonances, 
which require a full many-body theory for resonantly in- 
teracting Fermi atoms (li, [Ifl, [U, [23, [13, [H]. For a 
deep lattice this has been studied at zero temperature by 
Koetsier et al. [25| . 

We believe that the optical lattice introduces consid- 
erable new physics [3, [2^, [23, [H, [2^ , and our ambition 
here is to understand the resulting changes in the thermo- 
dynamics by considering a simple model for the system: 
the case of ideal particles (i.e. we neglect interactions 
for both the fermionic atoms and bosonic dimers) in a 
translationally invariant lattice potential. In this model 
the only effect of the Feshbach resonance is to maintain 
chemical equilibrium between the two species, while a 
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major simplification, this captures the essential physics 
of atom-molecule conversions in experiments [13, [3l| , and 
can be shown to be exact in the limit where the resonance 
is infinitey narrow [s^l • We note that for a (finite) narrow 
Feshbach resonance the inclusion of interactions should 
only impact our results quantitatively. 

In this paper we look at the entropic behavior of the 
system with special focus on the low-temperature regime, 
since isentropic surfaces characterize what will happen to 
a system in an adiabatic process. Such adiabatic process- 
es are often used in experiment, and our particular case of 
interest here is to predict how slowly ramping E^es across 
the Feshbach resonance will effect the temperature, and 
hence the degeneracy of the system. By looking at the 
temperature variation along such curves we can identify 
in which regimes cooling will occur, and surprisingly this 
turns out to depend on the average number of atoms per 
site in the optical lattice. Lattice physics plays an impor- 
tant role in changing the thermal properties of the sys- 
tem, indeed similar studies to what we present here have 
been undertaken for pure ideal Bose and Fermi gases in 
the translationally invariant lattice [33, 33| and more re- 
cently for the case with additional harmonic confinement 
[ssl . [3^ . [37I . Interactions play an important role in deep 
lattices, and more recent work has examined the effect 
these interactions have on the thermal excitation gener- 
ated durin g th e preparation of bosonic Mott-insulating 
states (ssl. l39l|. and on the feasibility of achieving the 
fermionic Neel state [43|. Finite temperature mixtures of 
atomic Bose and Fermi gases (with no conversion mech- 
anism) in lattices have been studied [4l| in an attempt 
to explain recent experiments (i^ . 



II. FORMALISM 

We consider a dilute system of Fermi atoms of mass 
TOa in a simple cubic optical lattice potential with 
sites and subject to periodic boundary conditions. The 
total number of atoms, A^tot, are divided equally into two 
different internal states which we denote as |t) and ||), 
with populations and , respectively. We define the 
filling fraction rj to be the average number of each type 
of atom on each site 
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(1) 



By applying a tunable magnetic field we can make the 
two atomic species interact with a Feshbach resonance 
and control whether these particle exist as individual (un- 
bound) atoms or participate in bosonic dimers with mass 
rrim = 2ma. We assume that the atoms and diatomic 
molecules are in thermal and chemical equilibrium, and 
we define the mean number of atoms (Na) and molecular 
dimers (Nm) with the number conservation condition 



A. Thermodynamics 

The thermal equilibrium condition for our system en- 
sures that the atoms and the molecules share the temper- 
ature T, and that the chemical potentials of the atoms 
and molecules are related as fim = + • Since we are 
considering an equal spin-mixture, this simplifies to 



(3) 



where = ~ fi*^ is the common atomic chemical po- 
tential, and fim is the molecular chemical potential. As 
the atoms only interact through the Feshbach resonance, 
we can describe the system by the single-particle energy 
levels, and the occupation of these are given by the Fermi- 
Dirac and Bose-Einstein distributions for the atoms and 
the molecules, respectively. Taking into account that we 
can effectively shift the molecular energy spectrum by 
an amount E^cs relative to the atomic spectrum via the 
Feshbach resonance, the relevant distributions are 
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(4a) 
(4b) 



The quantities E^ and E^^^ denote the single particle ener- 
gy levels in the lattice for the atoms and molecules respec- 
tively, where r is an appropriate quantum number (r = 
0,1,2,...). Since the molecules are governed by Bose- 
Einstein statistics we have the usual constraint that the 
(molecular) chemical potential must lie below the low- 
est molecular single particle state, i.e. fim < E^ + i?rcs- 
For each choice of the resonance energy, filling fraction, 
temperature and lattice depth the chemical potential /j, is 
then determined from the conservation of the total par- 
ticle number ^ where 



r=0 
00 

^rii ~ ^ ] fm- 



(5a) 
(5b) 



(2) 



Once the relation between ^, T and the externally ad- 
justable parameters has been established, it is straight- 
forward to calculate any thermodynamic quantity such 
as the entropy [43| 

Stot ^ Sa + Sm, (6a) 

00 

Sa=^^2j2 U: In /: + (1 - ,/:) ln(l - /^l fe, (6b) 

oo 

= - E [/™ /™ - (1 + 1^(1 + frn)] , (6c) 

where S'tot is the total system entropy, which is addi- 
tively formed from the atomic (Sa) and molecular (Sm) 
subsystem entropies. 
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B. Energy levels in the lattice 

The lasers which create the optical lattice have the 
wavelength \l = ^^t^I^l and the resulting potential 
is V^(x) = Vo,o-i(x) with a = a,m, for atoms and 
molecules, respectively, and where 

L(x) = sin^(fcLa;) + ii\Ti'{ki,y) + svTi'iki^z) , 

is the dimensionless shape of the lattice. The potential 
depth Vo^cr for the two different species is not the same, 
since a molecule consists of two atoms and therefore ex- 
periences twice the Stark shift of an individual atom: 
Vo.m = 2Vb,a = 2Vb. 

To find the possible energy levels we solve the station- 
ary Schrodinger equation, which for our simple cubic lat- 
tice separates into three ID Bloch equations. To find the 
3D energies it is therefore sufficient to calculate the ID 
energy levels and make all possible triplets thereof. It 
is convenient to rescale the stationary Schrodinger equa- 
tion for the species a by the corresponding recoil energy 
En^, = n^kl/2my. 

+ K),.i(x)| VA^) = Kri'U^), (7) 

where the barred quantities are energies in units of Eji^a- 
The operator on the left hand side only depends on the a 
through the species-dependent lattice potential Vb^o- and 
therefore the spectrum can be written as = E^CVo^a) 
for some species-independent function E^ . Furthermore, 
since the molecules are twice as heavy as the atoms and 
feel twice as deep a potential, we have the relation Vo,m = 
4Vo,a, which means that the energy spectra for Vb,m and 
4Vb_a have precisely the same shape but with different 
scalings on the energy axis. 




Figure 2: The density of states for a 31 x 31 x 31 optical lattice 
as a function of the potential depth for the particle type a. 
The density of states is found by binning the energy levels in 
intervals of width O-OAEr^^. 



1. Energy scales 

The fact that the atoms and molecules experience dif- 
ferent lattice potentials thus leads to a relative shift of the 
two spectra in addition to the magnetic field adjustable 
Feshbach detuning -Bres- Consequently, the position and 
width of the Feshbach resonance, as indicated by the in- 
terconversion of atoms and molecules, depends on the 
depth of the optical lattice potential [4j|. 

In Figure [2] we illustrate the energy spectrum in the 
3D lattice by plotting the density of states as a function 
of the lattice depth. The energy levels are distributed 
in a band-structure, where the bandwidths and the band 
gaps are strongly dependent on the depth of the lattice. 
The width of the allowed energy bands decreases as the 
lattice becomes deeper and for very deep lattice poten- 
tials the lowest energy band approximately reduces to a 
single energy level with a degeneracy that scales with the 
size of the system. On the other hand, the gaps between 
energy bands increase as the potential depth increases. 

Because of the different potentials the atoms and 
molecule experience, there is an energy splitting between 
the lowest atomic and molecular levels in the optical lat- 
tice. Using a tight binding analysis, appropriate to the 
regime Vb > AEji^a, we find that this splitting is given by 

as derived in Ref. (43 |. 

In a ID periodic potential a band gap opens at the 
edges of the first Brillouin zone for an arbitrarily small 
lattice depth. On the contrary, in a 3D system the contin- 
uum is only broken up into separated bands if the lattice 
potential is sufficiently deep, as is clear from Figure[2l In 
general it is of interest to know the width of the ground 
energy band {S^ and the gap between the ground and 
first excited energy bands {A^ in the lattice. We find 
that 

where both results are valid in the tight binding regime 
with the first result derived in Appendix El and the sec- 
ond result is obtained from the harmonic oscillator ap- 
proximation for the tunneling matrix element. We re- 
mark that Aa-/kB sets the temperature scale for the ex- 
cited bands to be relevant for the thermodynamics of 
particle to type a in the system. 

A convenient energy scale is the Fermi energy 
Ep(r], Vb), which is taken to be the highest occupied en- 
ergy level, when all the atoms are unbound. We note 
that for our choice of energy origin the relevant Fermi 
temperature for characterizing degeneracy is given by 
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Figure 3: The system entropy S'tot in the {ErcB,T) plane for 
Vo = WER,a and 77 = 0.8. 



Tp = (E-p — E^)/kB, with the atomic ground state 
energy. If 77 < 1 ther Fermi energy lies in the ground 
band (for sufficiently low temperatures) and the degen- 
eracy condition is fceT < Sa- 



III. ENTROPY 

To characterize the state properties of our system we 
numerically calculate the chemical potential for various 
points in the phase space (Vq, ?7, i?rcs, as discussed in 
Sec. ini Of most interest to us here is the behavior of 
the system entropy ^ as E^-qs and T vary, for which 
a typical example is shown in Figure [H We make the 
following general observations on this dependence: 

Atomic regime: Well above the resonance, i.e. for 



Er, 



> 



2Ep, the entropy is independent of E^ 



In this regime the lowest molecular energy levels 
are at a much higher energy than the atomic lev- 
els and are inaccessible. This result is only valid in 
the regime where ksT < Ep, so that we can neglect 
thermal activation of the molecular states. 

Molecular regime: Well below the resonance, i.e. for 
Eics ^ —2Ep, the opposite regime is entered in- 
to, in which only the molecular states are accessi- 
ble at low temperatures. Hence, in the lower left 
hand corner of the phase diagram the entropy is al- 
most independent of E^cs ■ At higher temperatures 
the resonance energy plays a larger role, since the 
molecules start dissociating at a temperature on 
the order of |i?ros|/fcB. 

Entropy peak: A notable feature of the entropy phase 
diagram is the behavior at E^^s ~ ^-f where the 
entropy as a function of the resonance energy has 
a maximum along each isotherm - a feature most 



noticeable at low temperatures. We investigate 
this feature further below, but note that it aris- 
es where the system is an equal mixture of atoms 
and molecules, and thus able to obtain maximum 
disorder for a given temperature. 



A. Entropy peak 

To analyze this last effect in more detail we look at 
the contributions to the entropy from both the atoms 
and the molecules along lines of constant temperature for 
^bT < 0.1i?F. We will show that when approaching the 
transition zone, where the conversion between atoms and 
molecules takes place, the entropy increases from either 
side. 

To clarify the underlying physics we develop a simple 
analytic model that should provide a good description in 
the deep lattice limit where the lowest energy bands are 
fiat, provided the temperature is small compared to the 
first band gap but large compared to the ground band- 
width, i.e. 



(11) 



We also require that ?7 < 1 so that higher bands are not 
populated by atoms due to the Pauli exclusion principle, 
and furthermore that a condensate does not occur in the 
molecular system, which would require a unique ground 
state. We note that these conditions are broadly consis- 
tent with the typical regime that experiments operate in. 
In this case we can approximate the density of states for 
the atoms and molecules as 



Pa{E)^2M'S{E~E',), 
p,n{E)^M^S{E-El-Eros), 



(12a) 
(12b) 



where, in the validity regime of this model, we can use 
([H]) to relate the difference between E^ and i?^^ to the 
lattice depth. 

In this approximation the level occupations (of the 
ground band) are independent of the particular level un- 
der consideration, and depend only on whether the par- 
ticle is an atom or a molecule, for which we denote the 
mean level occupation as /° and respectively. The 
corresponding entropy contributions will be 

Sa = -2M^ [f° In /o + (1 - /O) ln(l - /«)] fee, (13a) 
-M' [/° In/^ - (1 + ,C) ln(l + fj] fcB. (13b) 

We prove in Appendix [B] that at fixed temperature the 
total entropy S'tot = Sa+S„i, combining results (|13ap and 
(|13bp , is maximized when the resonance energy takes the 
value i?res""' — "^^a ~ ^^"i ill the deep lattice limit 
this may be approximated by 
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Figure 4: The atomic and molecular contributions to the entropy as well as the total entropy as functions of iSrcs at constant 
temperature, fceT = 0.04£f, lattice depth Vo = 10i?fl,a and for two different filling fractions. The horizontal dashed lines 
are the atomic and molecular entropy plateaus, Eqs. I|15ap and l|15bp . respectively, and the vertical dashed lines indicate the 
resonance energy maximizing the entropy as approximated by l|14p . 



using ([HI) and setting i?" ~ E-p. This coincides with 
the resonance energy at which the system obtains a 50% 
molecule fraction at zero temperature [43 |. 

In Figure [4] we present two examples from the numer- 
ical analysis of the behavior of the system entropy along 
an isotherm, and the portion of the system entropy aris- 
ing from the atomic and molecular subsystems. In both 
cases it is clear that the system entropy is maximized in 
the transition region at a value of the resonance energy 
that is in good agreement with our analytic expression 
for iSre""". The first case of 77 = 0.5 shown in Figure H] re- 
veals that the maximum in entropy arises from the rapid 
growth of degrees of freedom as population is coupled 
into the new subsystem (i.e. the growth of the atomic 
(molecular) entropy as we approach the transition from 
the left (right)) exceeding the reduction of entropy from 
the subsystem loosing population. 

In the second case of 77 = 0.9 shown in Figure [U we 
see behavior not fitting the previous description: not on- 
ly the molecular, but also the atomic subsystem entropy 
increases as the transition is approached from the right. 
The key feature leading to this is the non-monotonic de- 
pendence of the atomic subsystem entropy on ii^rcs, which 
occurs when the filling satisfies t] > \. We note from 
(fT3|) that Sa vanishes for /° = and = 1, and is 
maximum for = \- Thus for high filling (?7 > ^) as 
the resonance energy is lowered and the transition region 
is approached (from above) /° decreases, initially lead- 
ing to an increase in S'a, until /° decreases below i, at 
which point Sa begins to decrease. In contrast, due to 
their bosonic statistics the molecules give a contribution 
to the total entropy, which increases with increasing /,° . 

This peak in the atomic entropy can also be understood 
from l|B4a| . which states that the sign of dSa/ dE^^^ is 
determined by the sign of E^ — /j,, and this can be nega- 
tive if the chemical potential lies higher than the ground 
atomic energy level. This requires ^ < rj < 1, since in the 



deep lattice limit fi = E^ — /cbTIh 



V 



on the atomic 



side of the resonance, assuming ifTTj) . Note that there is 
no analogue effect for the molecules, since we always have 
Em +-E'res < 2/i in (|B4bp due to Bose-Einstein statistics. 



B. Entropy plateaus 

On the atomic side of the transistion zone {E^es ^ ZSp) 
and for the range of T where (fTTj) is valid, fa is constant 
at the value of rj. The associated entropy is entirely due 
to the atomic subsystem and is constant at the plateau 
value 



J,, 



'2M^ [?/lnr;+ (1 - 77) ln(l - t])] /cb. 



(15a) 



Also, in a similar temperature regime on the molecular 
side of the transistion zone (ii^ros S — ZSp) we have a 
constant molecular occupation /,„ = 77 and hence an as- 
sociated entropy plateau of 



-M^ [77 In ?/ - (1 + 77) ln(l + 77)] , 



(15b) 



These plateaus are indicated in Figure |4] as E^cs varies 
and clearly provide a good description of the total sys- 
tem entropy on either side of the transition region. In 
Figure [5] we show the total entropy for the system as a 
function of temperature for i?res = +4:Ef (in the atom- 
ic regime) and E^os = — 4i?F (in the molecular regime), 
respectively. These results clearly show the importance 
of the energy scales 5a and Ao-: For fceT ^ 

the plateau is observed in good agreement with the an- 
alytical prediction. For lower temperatures, k^T ^ 6^, 
the ground band is not uniformly occupied leading to a 
sharp suppression of entropy near T = 0. At higher tem- 
peratures. Act ^ k^T, excited bands become thermally 
accessible and contribute additional entropy. 
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Figure 5: Entropic behavior for Vo = 10Jf?fl,a. Left: The total entropy at i?ros = —4Ef (red curve) and Eics = +4i?F (blue 
curve), corresponding to a pure Bose and Fermi gas, respectively. The curves have plateaus, which are marked by the two 
dashed horizontal lines. Right: The height of the plateaus are read off manually (□ for the molecular and o for the atomic 
entropy plateau). These are well approximated by Eqs. I|15ap and ijlSbp (dashed curves) and intersect at the filling fraction 
rjc = 0.726. 



As is clear from the results in Figure [5] in the tempera- 
ture range So- <C ^ the entropy on the plateaus 
can be quite accurately determined from the full numer- 
ical calculations. We compare the results of such values 
against the analytic expressions (fT5|) for a range of filling 
fractions in Figure [5] on the right. The curves for the 
atomic and the molecular entropy plateaus intersect at 
the filling fraction rjc 

Tic = 0.726, (16) 
as is seen by equating (|15ap and (|15bp . 
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Figure 6: The dependence of the final temperature on the 
initial temperature in adiabatic Srcs-sweeps from the atomic 
to the molecular regime for different filling fractions and for 
a lattice depth Vo = 20iJfl,a. For r} < rjc the temperature 
increases while it decreases (for sufficiently low Ti) for rj > rjc. 



IV. ADIABATIC SWEEP OF E, 



Because ultra-cold atom systems are effectively iso- 
lated, i.e. not thermally connected to a reservoir, the 
various manipulations that can be made to the sys- 
tem will usually result in a change in its equihbri- 
um temperature. Such changes in temperature that 
occur during loading into optical lattices can be ap- 
preciable and has been of considerable recent interest 
[33, S, M, \M, [33, m, [li, US, SH- Here we investi- 
gate the temperature changes arising in the Feshbach 
coupled Fermi gas in an optical lattice. To do this we 
consider our initial system in the atomic regime with 
temperature Ti and entropy Si(Ti) and then determine 
the final temperature after an adiabatic sweep of i?ros 
into the molecular regime by equating the entropy func- 
tions, Sf{Tf) = Si{Ti), where Sf is the system entropy 
at the final i?ros value. Figure [6] shows our results for the 
relationship between Ti and Tf for various filling frac- 
tions. For moderate fillings, rj = 0.4,0.6,0.7, we observe 
that the sweep leads to a final temperature much greater 
than the initial temperature. This can be qualitatively 
understood because there are half as many molecules as 
atoms, and their temperature needs to be higher to carry 
the same amount of disorder as the unbound atoms. Of 
course this argument is not general. As discussed ear- 
lier, in the atomic regime in conditions where only the 
ground band is accessible, as the filling fraction increases 
past 7] ~ 0.5, the entropy of the fermionic atoms decreas- 
es (at fixed T) - on the other hand, the bosonic entropy 
only increases with increasing rj. Our results show that 
for the parameters under consideration in Figure [6] that 
for 77 > 0.8, there is a wide region of initial temperatures 
for which the Feshbach sweep leads to a much colder 
molecular gas. 
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Figure 7: Schematic illustration of the adiabatic cooling. Left: If the atomic entropy plateau lies above the molecular one 
iv < ')c) the final temperature in an adiabatic -Eros-sweep will be higher than the initial temperature as shown to the left. 
However, if the molecular entropy plateau lies above the atomic one {rjc < i] < 1) the final temperature will be comparable to 
the bandwidth of the molecular spectrum, Sm- 



This behavior is quantitatively explained in Figure [7| 
where the entropy curves in the atomic and molecular 
regimes are shown for two different cases: (1) If the atom- 
ic entropy plateau lies above the molecular one {rj < rjc) 
and we start with an initial temperature Ti correspond- 
ing to an atomic entropy on the plateau, then Tf will 
be independent of and lie above Ti. (2) If on the other 
hand the atomic entropy plateau lies below the molecu- 
lar one {r/c < rj < 1), the final temperature will be much 
lower, i.e. lie in the regime where fceT is comparable 
to Sa- A detailed prediction for the final temperature 
requires a theory for the molecular entropy dependence 
on temperature outside of the plateau region. There is 
no general closed-form expression for this in the lattice 
potential, however in the regime where significant tem- 
perature reduction is predicted the uniform Bose gas ex- 
pression should be applicable with the replacement of 
mm by an effective mass including the infiuence of the 
lattice (it should be noted that in the molecular regime 
for kBT < Sm condensation is likely and will need to be 
included to make a quantitative prediction). 



The above results only apply when the the initial 
state lies on the entropy plateaus (i.e. E^cs ^ 2i?F and 
Sa k^T <^ Aq). The upper temperature limit, Aa, 
increases with increasing Vb, and the lower limit, 5a, de- 
creases with increasing Vb, so that the region of appli- 
cability of this simple analysis improves with increasing 
lattice depth. However, this is also the regime in which 
interaction effects become more important requiring a de- 
scription beyond that presented here, if quantative pre- 
dictions are to be obtained. 



V. THE ROLE OF INTERACTIONS 

We have thus far neglected interactions in our analysis 
to elucidate the aspects of the thermodynamics arising 
from the band structure in the lattice potential. We now 
briefly discuss how interactions will modify our conclu- 
sions. Our analysis of isentropic sweeps of the resonance 
energy relied on the existence of plateaus of constant en- 
tropy on both sides of the Feshbach resonance. We there- 
fore concentrate on how interactions impact the entropy 
in the molecular and the atomic limits. 

On the molecular side of the resonance we have a pure- 
ly bosonic system with weakly repulsive residual interac- 
tions. As the lattice depth is increased correlation effects 
get stronger in the system. For the translationally invari- 
ant lattice the special case where the number of bosons 
is commensurate with the number of lattice sites then 
the system can undergo a quantum phase transition to 
the Mott-insulating phase. In the Mott-insulating phase 
a gap, given by the on-site interaction strength (J7m), 
emerges in the excitation spectrum. For temperatures 
larger than Um/ks (but much smaller than the gap to 
the flrst excited lattice band ) th e molecular entropy will 
still exhibit a plateau (45l.[4^.l4^. As in the ideal case the 
existence of an entropy plateau is a consequence of the 
uniform distribution of the particles over all the states 
in the lowest band. Hence the value of the molecular 
entropy on the plateau is determined mostly from com- 
binatorial arguments. At lower temperatures the exci- 
tation gap exponentially suppresses the entropy of the 
molecules. 

It is worth noting that the main results of this work 
have been in the regime rj < 1, for which the boson num- 
ber is incommensurate with the number of lattice sites 
and the Mott-insulator transition can not occur [i^ . For 
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this case the our ideal calculations should have a wider 
regime of validity. 

On the other side of the resonance the interaction be- 
tween the free atoms is attractive and the ground state 
is a paired BCS superfluid, characterized by a pairing 
gap, Apair. The size of the pairing gap depends on 
the strength of the attractive interaction Ua compared 
with the bandwidth Sa- In the weak coupling limit Apah- 
is exponentially small and proportional to Sa, while for 
stronger interactions the size of Apair is set by \Ua\ [49| . 
As was the case for the molecules the atomic entropy is 
reduced by an exponential factor for smaller than 
Apair, while it reaches a constant plateau in the inter- 
mediate temperature regime Apair < ^ Aq. Hence 
correlated states introduce a new low energy scale |J7o-|, 
which replaces the bandwidth Sa as the relevant energy 
scale for the existence of an entropy plateau, and our con- 
clusions remain qualitatively correct when interactions 
are included, provided \Ua\ <C A^ and U,n ^ A,„ in the 
deep lattice limit. Whether this inequality can be satis- 
fied will depend on the details of the interactions. 

Finally, we remark that while the Feshbach molecules 
are produced in the least bound ro- vibrational state they 
are remarkably robust against collisional deexcitation, 
which could lead to losses when a molecule shares a lat- 
tice site with another molecule or a free atom. This sta- 
bility arises from Pauli blocking of the the constituent 
fermionic atoms [HOl- Again we note that since our pri- 
mary interest here is in systems with fiUing fractions less 
than unity, the effect of coUisional losses will be smaller. 

When an external harmonic trapping potential is 
present the qualitative similarity of the physics of the 
ideal and strongly correlated regimes is not clear. For in- 
stance, in this case the gas will in general contain several 
distinct spatial re gion s exhibiting either superfiuid or in- 
sulating behavior [5l|. Hence this situation will require 
future work. 



VI. CONCLUSIONS 

In this paper we have considered the Feshbach reso- 
nance based association of a Fermi gas of atoms con- 
fined in a uniform optical lattice potential. Our anal- 
ysis of the statistical mechanics of this system has re- 
vealed a rich range of behavior in this system arising from 
the strikingly different degenerate properties of the atom 
and molecular degrees of freedom, and the well-isolated 
ground band that forms in moderately deep lattice poten- 
tials. We have considered the effect of an adiabatic ramp 
of the Feshbach resonance from the atomic to molecu- 
lar regimes, and shown that a colder molecular gas can 
be produced for a wide parameter regime if the filling 
fraction lies in the range 0.726 < jy < 1. 

By focussing on the ideal gas case we have been able 
to clarify the effect of the translationally invariant lattice 
potential on the atom-molecule equilibrium. Building on 
this study, future work will be to include interactions and 



the additional confinement present in experiments due to 
the inhomogeneous harmonic trap potential. 
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Appendix A: ANALYTIC EXPRESSIONS FOR 
THE FIRST BAND GAP 

In the deep lattice limit we can make a tight-binding 
approximation and regard the potential as a collection of 
independent wells. The central band energy is then 
given by the discrete energy levels E'^ in each well which 
to lowest order are harmonic oscillators, but anharmonic 
corrections can be calculated using perturbation theory 
(see Appendix A of Ref. Coupling between wells 

induces a finite bandwidth, that we neglect in the deep 
lattice (fiat band) limit. Because the lattice potential is 
separable the energies are of the form 



(Al) 



where we refer to the as the ID energies along the 
Xj direction, and {n} {nxjUyjUz} are the quantum 
numbers describing the single well state. The lowest two 
ID energy levels are given by [43 | 




1 



Vo,a 5 49 1 




From (jAip we see that the lowest 3D energy band lies at 
= while the next 3D energy band lies at E^ = 
2£^ + £^ yielding a bandgap of 



(A3) 



A, = Ei - E''= £l - £° 



Approximating the two lowest ID energy levels by (jA2p 
we have 




1 - 



23 



1 



24 



Eb.,.- (A4) 



Appendix B: MAXIMISING ENTROPY IN THE 
FLAT-BAND LIMIT 



We are interested in how the total entropy Stot = Sa + 
Sm changes if we change the resonance energy but keep 
the temperature constant, i.e. 



f dSt, 



\dE^., 



dSa ( dfa 



■cs / T 



res / T 



dSm ( df^^ 



(Bl) 



■OS / T 
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The number of unbound atoms is Na = 2M^/° and the 
number of molecules is Nm = M^f^, so by combining 
([T]) with the condition of particle conservation ^ we get 



fa + fl - V- 

We may thus rewrite I|B1[) as 



dE, 



res / T 



(B2) 



(B3) 



I CS / T 



From (|T3|) the partial derivatives of Sa and S"™ with re- 
spect to the occupation numbers and can be found, 
leading to 



dSm ^ M''{El+E,,,~2^i) 

a/0 T 

which we can insert in (|B3p . yielding 
/ d5tot \ _ {2El -El- i?,es) ( df^ 



dE, 



res / T 



(B4a) 
(B4b) 

(B5) 



The dependence of E^cs on /o is via the chemical po- 
tential, so the derivative is 



dfl 

dE,.cs / T 



'AE°-^i)/kBT 



(B6) 



which is positive since d^/dE^^^ > in this regime |44 |. 
Thus we see from (|B5P that the sign of the derivative of 
S'tot with respect to £^res is the same as the sign of the 
quantity 2E^ — E^ — i^res- We conclude that the entropy 
is maximal when the resonance energy equals 



E^-- ^ 2El -El^ El - {El - El) 



(B7) 



and in this deep lattice limit we can approximate E^ — El 
by ([H]) and set E^ equal to the Fermi energy, such that 
we get 



E" 



Ev- -- 



1 



16 ^Vo,a/En,a 



E 



R,a- 



(B8) 



res / T 
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